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ABSTRACT 

Within the magnetar scenario, the "twisted magnetosphere" model appears very 
promising in explaining the persistent X-ray emission from the Soft Gamma Repeaters 
and the Anomalous X-ray Pulsars (SGRs and AXPs). In the first two papers of the se- 
ries, we have presented a 3D Monte Carlo code for solving radiation transport as soft, 
thermal photons emitted by the star surface are resonantly upscattered by the magne- 
tospheric particles. A spectral model archive has been generated and implemented in 
XSPEC. Here we report on the systematic application of our spectral model to different 
XMM~Newton and INTEGRAL observations of SGRs and AXPs. We find that the 
synthetic spectra provide a very good fit to the data for the nearly all the source (and 
source states) we have analyzed. 

Key words: Radiation mechanisms: non-thermal - stars: neutron - X-rays: stars. 



o 

O 



X 



1 INTRODUCTION 

Over the last few years, increasing observational evidence 
has gathered in favour of the existence of "magnetars" , i.e. 
neutron stars (NSs) endowed with an ultra-strong magnetic 
field {B « 10" - lO^'^ G), about ten times higher than 
the critical threshold at which quantum electro-dynamical 



(QED) effects become important {Be 



4.4 X 10" G). 



The family of magnetars candidates comprises two classes of 
sources, the anomalous X-ray pulsars (AXPs) and the soft 7- 
ray repeaters (SGRs). About fifteen objects are known, all 
characterized by similar properties: slow X-ray pulsations 
(P ~ 2-12 s), large spin-down rates (P ~ 10""-10"^" s/s), 
atypical persistent X-ray luminosity of ~ lO'^^-lO'^® ergs~^, 
lack of bright optical companions (favouring an interpreta- 
tion in terms of isolated objects), and a high level of burst- 
ing/ flaring activity which can differ among the two classes 
fsee lWoods fc Thompsonll2006l : lMereghettill2008l . for recent 
reviews.) 

Spectral data can provide key information about the 
physics of ultra-magnetized NSs. High-energy observations 
of SGRs and AXPs persistent emission cover now the ~ 0.5- 
10 (XMM-Newton / Chandra) and ~ 20-200 (INTEGRAL) 
keV bands. It is possible, although not proved yet, that dif- 
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ferent emission mechanisms are responsible for the emission 
in the two bands. 

The low-energy (< 10 keV) spectrum of the persistent 
(i.e. outside bursts) emission is often modelled in terms of a 
blackbody (kT ~ 0.3-0.6 keV) plus a power-law with photon 
index Fs ~ 2-4, although in some AXPs a two blackbody 
model has been also applied. The X-ray persistent emis- 
sion above 20 keV has a power-law spectral shape (Fh ~ 2) 
which, in particular in AXPs, is markedly harder t han that 
observed below 10 keV (see again 'Mereghettil l2008l . and ref- 
erences therein). However, although these phenomenological 
fitting models have been systematically applied over the last 
decade, a convincing physical interpretation of the various 
spectral components is still missing. 

In the magnetar framework, the t wisted magnetosphere 
mode l ( Duncan fc ThompsonI 19921: Thompson fc DuncanI 
1 19931 : [Thompson. Lvutikov fc Kulkarnil 120021 ') offers a 
promising physical interpretation for the X-ray emission 
from SGRs/AXPs. In particular, the supporting currents 
required to sustain the twist are substantially in excess of 
the Goldreich- Julian current and can produce large optical 
depth to resonant cyclotron scattering (RCS). Soft (ther- 
mal) photons produced at the star surface gain energy in 
repeated scatterings and this leads to the formation of an 
extended, high-energy tail, superimposed to the seed ther- 
mal component. The qualitative predictions of this model 
have been verified to match some spectral and timing prop- 
erties of magnetar sources, as the spectral shape observed 
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in quiescence below ~ 10 keV a nd the long term varia- 
tion observed in some sources (e.g. Mereghetti et al.ll2005bl : 
iRea et al.ll2005l : ICampana et al.l[2007f K 

More recently, several efforts have been carried out re- 
cently in order to test the model quantitatively against real 
data using different approaches and with a varying degree 
of sophisticat ion. A first attempt in thi s direction has been 
presented by iLvutikov fc Gavriill (|2006t ). These authors de- 
veloped a very simple, semi-analytical treatment of the RCS 
process by working in 1-dimensional geometry. They as- 
sumed that seed photons are emitted by the NS surface with 
a blackbody spectrum, and Thomson scattering occurs in 
a thin, plane parallel magnetospheric slab permeated by a 
static, non-relativistic, warm medium at constant electron 
de nsity. These mode ls have been implemented in XSPEC 
bv lRea et ahl (l2007al lbl. [200i 'l and successfully applied to all 
magnetars spectra below lOkeV. The good agreement found, 
even within a simplified treatment, supports the idea that 
RCS in a sheared magnetosphere plays a central role is the 
formation of magnetar spectra in the (soft) X-ray range . 
The same RCS model has been used bv lGiiver. et al.l (|2007l ). 
who assumed that seed photons comes from an atmosphere 
surrounding the star. More recent ly, 3-D Monte Carlo cal- 
culati ons have been presented by [Fernandez fc ThompsonI 
l|2007l ). although these spectra have never been applied to 
fit X-ray observations. 

Motivated by this, we present here a systematic ap- 
plication of our 3D Monte C arlo spectral calculations (see 
iNobili. TuroUa fc Zanj|2008al lbl for aU details) to AXPs and 
SGRs spectral data. The paper is organized as follows. In 
§ [2] we briefiy summarize the basic features of the model. 
The data sample and spectral results are presented in § |3l 
Discussion and conclusions follow in §5[4l[5l 



2 THE MODEL 



As discussed in detail in iNobili. Turolla fc Zan i (|2008al lbl). 
we have recently developed a 3-dimensional treatment of 
RCS, aimed at a detailed investigation of the spectral, 
timing and polari zation properties of magnetars (see also 
iPavan et al.l 120091 ). Our Monte Carlo code, which is com- 
pletely general and can handle different magnetic field 
topologies and distributions of seed photons, has been used 
to produce an archive of spectral models that have been 
subsequently implemented in XSPEC. Such models rely on 
a number of choices fo r the assumed configura ti on, tha t are 
discussed in detail in INobili. Turolla fc Zand (|2008al ) and 
briefiy summarized in this section. 

First, the XSPEC archive models have been computed 
by assuming that the star surface emits as a blackbody at 
an uniform temperature, kT, and that the surface radia- 
tion is isotropic and unpolarized. The magnetic field topol- 
ogy is assumed to be a twisted, force free dipole and is 
uni quely characterized by the value of the twist angle, A0 
(see [Thompson. Lvutikov fc Kulkarnill2002l ). No attempt is 
made to fit the value of the polar field strength, that has 
been fixed at 10^* G. Our model is based on a simplified 
treatment of the charge carriers velocity distribution which 
accounts for the particle collective motion, in addition to 
the thermal one. Pair production has been neglected. Mag- 
netospheric electrons stream freely along the field lines (the 



motion across the field lines is quantized). The electron ve- 
locity distribution parallel to the field is taken to be a 1-D 
relativistic Maxwellian at temperature Te, superimposed to 
a bulk motion with velocity Pbuik (in units of the light ve- 
locity c). Besides, in order to minimize the number of free 
parameters, the models in the archive were computed as- 
sum ing that the electron temper ature is related to Pbuik 
(see lNobih. Turolla fc Zanel l2008a. for aU details). Scatter- 
ing in a magnetized medium was treated by considering only 
the resonant part of the magnetic Thomson cross section 
and neglecting electron recoil along the field direction. For 
the sake of conciseness, in the following this approximation 
will be referred to as the (resonant) Thomson limit. On the 
other hand, the code is completely general and inclusion 
of the relativistic QED resonant cross section, which is re- 
quired in the modelling of the hard (~ 20-200 keV) spectral 
tails observed in the magnetar candidates, i s under way (see 
INobili. TuroUa fc Zanj l2008bl : jPavan et all '2009 ) . Electron 
recoil in the direction parallel to the field starts to be im- 
portant when the photon energy in the electron rest frame 
becomes comparable to the electron rest energy. If 7 is the 
mean electron Lorentz factor, this occurs at typical ener- 
gies ~ meC^/7. Assuming mildly relativistic particles, the 
previous limit implies that conservative scattering should 
provide good accuracy up to photon energies of some tens 
of keV. This, together with the fact that resonant scattering 
occurs in regions where B <C Bcru, makes the use of the 
(much simpler) non-relativistic (Thomson) magnetic cross 
section adequate. The final XSPEC atable spectral model 
(22 MB in size, named ntznoang.mod, hereafter NTZ) has 
been created by using the routine wf tbmd, available on- hneQ 
In summary, it depends on four free parameters {Ptuik, A</), 
log kT plus a normalization constant), which can be simulta- 
neously varied during the spectral fitting following the stan- 
dard minimization technique. In NTZ models the photon 
number is conserved and the monochromatic number fiux 
which reaches infinity is the same as that of the seed black- 
body spectrurr0- This implies that the normalization con- 
stant (norm) divided by kT^ is proportional to the emitting 
area on the star surface (see ^4.2|) . It is important to note 
that the NTZ model has the same number of free parame- 
ters than the canonical blackbody plus power-la w empirical 
mode l or the ID RCS model recently discussed in lRea et al.l 
(|2008f), and hence has the same statistical significance. In 
the following sections we present its systematic application 
to the soft (0.5-10 keV) X-ray spectra of magnetar candi- 
dates. 



^ see'http: / /heasarc.gsfc.nasa.gov/docs/he asarc/ofwg/docs/general/| 
modelfilGS_memo/modolfilos_memo.html. 

^ Photon number is not necessary conserved when Landau- 
Raman scattering i s ac counted for (photon spawning; 
INobili. Turolla fc Zanj [20P8b). Present models have been 
obtained under the assumption of conservative magnetic scatter- 
ing, so photon number conservation is ensured. Also, spectra in 
the model archive have been computed averaging over all viewing 
directions, i.e. photons are collected over the entire observer's 
sky. When one particular line-of-sight is chosen, the photons 
reaching the observer are clearly less than those emitted by the 
surface. 
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3 APPLICATION TO MAGNETAR'S SOFT 
X-RAY SPECTRA 

We applied the NTZ model to a large sample of AXPs 
and SGRs, using XMM-Newton and INTEGRAL data. 
T he sample o f data sets basically coincides with that used 
in iRea et al] (|2008| ). and we refer to that paper for more 
details and for a discussion on the data analysis. However, 
at variance with iRea et al.l ('2008), we did not consider here 
the transient magnetars XTE J1810-197 and CXOUJ1647- 
4552 , the analysis of which will be reported in a sep- 
arate paper devoted to a detailed investigation of their 
outburst evolution over a period of several years. On the 
other hand, we included the recent XMM-Newton obser- 
vations of CXOU JOlOO-7211 and SGR 1627-41 , that were 
not available at the time of the previou s investigation (see 
iTiengo et alll2008l : lEsposito et"aLl |2008| , for further details 
on these observations). 

All fits have been performed using XSPEC version 11.3 
and 12.0, for a consistency check. A 2% systematic error was 
added to the data to partially account for uncertainties in 
instrumental calibrations. Only in the case of CXOU JOIOO- 
7211 , which has a very low absorbtion, the fit has been 
performed in the 0.5-10 keV energy range. For the other 
sources, that are highly absorbed, the 0.5-1 keV energy 
range (0.5-2 keV energy range in the case of SGR 1627-41 ) 
was excluded from our spectral fittinglf] We then checked 
that, for all our targets, the values of A'^^f derived fitting 
the 1-10 keV EPIC-pn spectra, are consistent (within 1 a) 
with those obtained using the 0.5-10 keV range relative to 
the same data set. We us ed the more updated solar abun- 
dances by Lodders ( 20Qi ) , instead of the older ones from 
lAnders fc Grevessd l|l989t ). As a consequence, the value of 
the absorption is, on average, slightly higher than that re- 
ported in the literature for the same model. This does not 
affect the other spectral parameters. 

Tableland Figs.[T][2]report our results in the 1-10 keV 
range for those sources that do not show a substantial spec- 
tral variation, in which cases only one dataset (the longest 
available) has been considered. Also, in Tables [51 E] and Fig.O 
we show the fits for sources that do exhibit spectral varia- 
tion, in which cases we considered a set of two-three obser- 
vations for each source, corresponding t o different spectral 
states. At variance with lRea et al.l l|2008l '). in these cases the 
different datasets have been fitted by imposing that the ab- 
sorption, Nh, is the same. As it can be seen from the tables, 
in most of the cases we found that a NTZ model alone suc- 
cessfully reproduces the soft X-ray part of the spectrum up 
to 10 keV, without the need of further components (see SJH 
for a discussion) . With reference to the sample considered in 
iRea et al] (|2008l ). the only two sources for which we do not 
find a satisfactory fit are IE 2259+586 and 4U 0142+614 , 
which are discussed separately below. 

There are a few AXPs and SGRs in our sam- 
ple (namely IRXS J1708-4009 , IE 1841-045 , 4U 0142+614 , 
IE 2259+586, SGR 1900+14 and SGR 1806-20) that are 
known to have a conspicuous emission in the hard X-ray 
band. For these objects we have repeated our modelling in- 



^ The emission in this energy range is in fact mostly affected by 
interstellar absorption. M oreover, this is the band where most of 
the calibration issues lay l lHaberl et al.|[2"00'4h . 



eluding also their INTEGRAL spectra in order to inves- 
tigate whether our model can reproduce part of the emis- 
sion when considering the whole SEEQ. A further power-law 
(PL) has been introduced in order to account for the non 
thermal hard X-ray component, although at a purely phe- 
nomenological level. A free constant was multiplied when 
using both XMM-Newton and INTEGRAL data to ac- 
count for inter-calibration uncertainties (the values of the 
constant was always differing by less than 10% with respect 
to XMM-Newton which was set to unity). Results are re- 
ported in Tableland Fig.|4l We find that, again with the ex- 
ception of IE 2259+586 and 4U 0142+614 , in all other cases 
a NTZ+PL spectral decomposition successfully models the 
1-200 keV emission. However, the fit converges with a hard 
X-ray component that gives a substantial contribution when 
extrapolated to the soft X-ray band and, consequently, the 
best fitting parameters of the NTZ model are substantially 
different from those we found fitting the 1-10 keV emission 
only. 

As the tables show, in some cases the fit converges with 
a value of A(j) which is too close to the upper bound of 
our model archive {Acp — 2), making impossible to com- 
pute the whole (closed) la contour level. Indeed, in the 
case of the combined XMM-Newton and INTEGRAL fits 
of IE 1841-045 and SGR 1900+14, the parameters of the 
NTZ model becomes less constrained than when fitting the 
XMM-Newton data alone, and we could only set an upper 
limit on A(j) (see Table|4|. Although we do not regard this 
as a particular problem, we should caveat that in some fits 
the value of the twist angle appears to be less constrained 
than those of the other model parameters. 



4 DISCUSSION 

In this paper we have applied a 3D Monte Carlo model of 
resonant scattering to the phase averaged spectra of an ex- 
te nsive set of magn etar sources. The sample is the same as 
in iRea et al] (|2008l ). with the exception of the two tran- 
sient magnetars XTEJ1810-197 and CXOU J1647-4552 , 
and with the inclusion of recent XMM-Newton observations 
of CXOU JOlOO- 7211 and SGR 1627-41 . The discussion of 
our findings is organized as follows. In tj4.1l we first con- 
centrate on the results emerging from the fit of the XMM- 
Newton data only, i.e. on the energy band below 10 keV. 
We report the discussion of our spectral results, a compar- 
ison with the an alogous findings from the ID RCS fits of 
iRea et al.1 (120081 ). and we discuss the cases of IE 2259+586 
and 4U 0142+614 , i.e. the only two sources for which a fit 
with the NTZ model is unsatisfactory. A search for correla- 
tions among the sources and the model parameters, in the 
same energy band, is reported in ?I I4.2I In M4.3I we then dis- 
cuss the NTZ model in the contest of the interpretation of 
the whole SED up to ~200 keV, basing this time on both 
XMM-Newton and INTEGRAL data. The main hmitations 
and caveats are then summariszed in 94.41 



* In the case of SGR 1806-20 there are no observations taken at 
similar epochs in the hard and soft X-ray bands. For this reason, 
and because of the high variability of this source, no attempt has 
been mad e to fi t its broadband distribution either here nor in 
iRea et al.l 1I2OO8I') . 
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4.1 The soft X-ray spectra: fits of the 
XMM-Newton data 

Our results on the fits of tlie XMM-Newton data are re- 
ported in Tab. [1] and Figs. I1I2I for sources witliout a sig- 
nificant spectral variation, and in Tables (2] |3] and Fig. [3] 
for sources that do exhibit spectral variation. When we re- 
strict to the 1-10 keV band, we found that the NTZ model 
successfully reproduces the soft X-ray part of the spec- 
trum of most of the sources (apart from IE 22594-586 and 
4U 0142-1-614 ), without the need of additional components. 
This represents a substantial improvement with respect to 
previous attempts to model magnetars quiescent emission 
in the same ene rgy band with a simpler ID RCS model 
l|Rea et al.l |2008| ). where it was found that in a few cases 
a PL component was required, in addition to the RCS one, 
to provide an acceptable fit to the data below 10 keV. In 
particular, this was the case of few AXPs (including e.g. 
IRXS J1708-4009 and IE 1841-045 ) and of SGR 1806-20 , aU 
them detected also above > 20 keV. In such cases, we found 
that the ID RCS component reproduces the spectrum only 
up to 5-8 keV. In order to match the data at the highest 
XMM-Newton energies, the contribution of a PL must be 
included. On the other hand, the slope of this additional PL 
is the same of that describing the INTEGRAL spectrum in 
the ~ 20-200 keV band. It is well possible that the mech- 
anism responsible for the hard X-ray emission provides a 
non negligible contribution in the soft X-ray range. How- 
ever, our findin g is that the resonant C ompton scattering 
model bv lNobili. Turolla fc Zari3 (|2008al ) correctly describes 
the data in the whole XMM-Newton energy range, i.e. up 
to 10 keV, for all the sources (and source states) reported 
in Tables I1I2I3I The application of the NTZ model to the 
> 20 keV emission from magnetars is discussed later on. 
It is worth noticing that the NTZ model has t wo free pa- 
ramet ers less than the ID RCS-I-PL model used in lRea et al.l 
()2008l ). Hence, besides self-consistency, it is also more robust 
on a statistical ground. 

In all cases we found that Nh , as derived from the NTZ 
model, is lower than (or consistent with) that inferred from 
the BB-I-PL fit, and consistent with what derived from fit- 
ting the single X-ray edges (Durant & van Kerkwijk 2006). 
The same w as also true for th e fits with the ID RCS model 
presented in iRea et al.l l|2008h . This is not surprising, since 
the power-law usually fitted to magnetar spectra in the soft 
X-ray range is well known to overestimate the column den- 
sityo The surface temperature we derived fitting the NTZ 
model is slightly higher than the corresponding ID RCS 
temperature and consistent with the BB temperature in the 
BB-I-PL model. On the other hand, a quantitative compar- 
ison between the values of Ptuik and A(j}, i.e. the parame- 
ters that describe the magnetospheric currents in the NTZ 
model, with the corresponding parameters of the ID RCS 
model {(3t and Tres) is more difficult, due to the different as- 
sumptions about the currents velocity distribution and the 
magnetic field topology. The ID RCS assumes a plane par- 
allel slab (i.e. photons can only propagate along the slab 
normal, either towards or away from the star). Moreover, 

^ This is because the fitting procedure tends to increase absorp- 
tion (i.e. Nh) to counter the steep rise of the power-law at low 
energies, which eventually diverges as E — > 0. 



magnetospheric charges are assumed to have a top-hat veloc- 
ity distribution centered at zero and extending up to a given 
temperature, ifcTe. To some extent, this scenario mimics a 
thermal, ID, motion (in which case kTe can be assimilated 
to a mean e~ energy, i.e. to the temperature of the ID elec- 
tron plasma). Since the e~ velocity distribution averages to 
zero, no bulk motion is accounted for: if we only consider the 
frequency shift due to the e~ motion then a photon has the 
same probability to undergo up or down scattering. Photon 
boosting by particle thermal motion in Thomson limit may 
still occur, because the spatial variation of the magnetic field 
is taken into account. For a photon propagating from high 
to low magnetic fields, multiple resonant cyclotron scatter- 
ing will, on average, up-scatter the transmitted radiation, 
giving rise to the formation of an hard tail. Since, at vari- 
ance with the ID RCS model, our code accounts for both 
bulk and thermal motion, one may expect to find a value 
of Pbuik systematically lower than the value of the velocity 
parameter /9t in the ID RCS fits. However, we find that this 
is not always the case and the relation between the values of 
Pbuik and Pt appears to be more complex. The NTZ model 
does not explicitly contain the optical depth as a parame- 
ter. However, A(fi and Pbuik can be related to the (average) 
scattering depth, as discussed in iRea et al.1 (|2008h . In par- 
ticular, the value of Ta«e can be read off from their figure 10, 
simply by dividing the value of the depth corresponding to 
a given A(j> by Pbuik- The values of Tave we derive (~ 0.1- 
1, as expected since the average number of scatterings for a 
typical photon is of this order) are systematically lower than 
those obtained for the ID RCS Tres ■ This means that the ID 
RCS model is intrinsically less efficient in up-scattering the 
seed photons. Its spectrum is, in fact, softer and it requires 
a larger depth (and hence a larger number of scatterings). 
Moreover, the hardest sources require an additional PL to 
match the observed XMM-Newton spectra. When observa- 
tions of the same source at different epochs are available, 
it is of interest to check if the best-fit values of the model 
parameters change, since this can reveal how the physical 
properties of the surface/magnetosphere evolve in time. In 
the case of IE 1048-5937 all the parameters are compatible 
with being constant within the errors, with the only excep- 
tion of the normalization which appears to increase follow- 
ing the fiux rise. The same holds for SGR 1806-20 where 
there may be also an indication for a decrease in the twist 
angle after the giant flare of December 27th 2004. Errors 
are however quite large and prevent any firm conclusion at 
this stage. The time behaviour of IE 1547.0-5408 is puz- 
zling, since Pbuik increases while A<j) decreases as the flux, 
and the model normalization, increases. In this source there 
is also a quite robust evidence that the surface temperature 
increased with the flux. 

As we mentioned before, in the case of IE 2259-1-586 
and 4U 0142+614 we were unable to find a satisfactory fit 
with the NTZ model. We noticed that, below 10 keV, these 
two sources are characterized by a rather soft spectrum: in 
the canonical BB-I-PL decomposition the power law index is 
~ 4 instead of ~ 2 — 3 as in other magnetars. At the same 
time, the power law tail starts very close to the energy at 
which the BB peaks. Such a spectral shape appears difficult 
to explain in terms of upscattering of soft photon, what- 
ever the nature of the comptonization process might be. In 
fact, a PL component that starts close to the BB peak is 
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a signature of a full fledged comptonization in which case, 
however, it is expected to be quite flat. Conversely, a steep 
power law tail is typical of weak comptonization and depar- 
tures from the BB spectrum occur at energies beyond the 
peak. One possible explanation is that the BB peak appears 
to be less prominent to the observer because the region that 
emits th e soft seed photons is not com ing completely into 
view. In iNobih. TuroUa fc Za"n3 i 2008al ) ■ we considered the 
effects of a non-homogeneous surface temperature distribu- 
tion, by examining the case in which photons are emitted 
by a single surface patch. In Fig. [5] we show the results of 
our simulations in the case in which the emitting region is 
confined to an equatorial strip (left panel) or to a polar cap 
(right panel). The subdivision of the star surface, that of 
the sky, the energy range and bin width are the same as 
that used in Nobili, TuroUa & Zane (2008a; see Fig 7 in that 
paper). The different curves show the emerging spectrum, 
as viewed by an observer whose line of sight (LOS) makes 
an angle Os with the spin axis and for different values of 
the observing longitude, = 20°, 140°, 220°. These three 
values correspond to having the emitting surface patch in 
full view (seen nearly face on), partially in view and almost 
screened by the star. As it can be seen, when the emitting 
patch is in full view the observed spectrum consists of a well 
visible thermal component and an extended power-law-like 
tail. On the other hand, if the emitting region is not di- 
rectly visible, no contribution from the primary blackbody 
photons is present: the spectrum, which is made up only 
by those photons which after scattering propagate "back- 
wards" , has a depressed thermal peak and a much more dis- 
tinct non-thermal shape. In particular, in the examples of 
Fig. [5l the scattering efficiency is not very high and curves 
corresponding to $s ~ 220° have a steep power-law tail 
with photon index ~ 4, as observed in IE 2259-1-586 and 
4U 0142-1-614 . Although a quantitative fit of the spectra of 
these two sources including different thermal maps in the 
code would be unfeasible (because it would introduce too 
many degrees of freedom), it is tempting to speculate that 
the peculiar spectrum of IE 2259-1-586 and 4U 0142+614 is 
due to a strong inhomogeneity in the surface temperature 
distribution, with the hotter region almost antipodal with 
respect to the observer. This is also compatible with the 
fact that these two sources have a rather low pulsed fraction 
with respect to other magnetars. On the other hand, the 
double peaked pulse profile of IE 2259-1-586 is difficult to 
explain in this picture. Another possibility is that the phase 
average spectrum appears to be quite soft because it re- 
fle cts the contribution o f diff erent components. As discussed 
bylmodseFaO llooi) and iRea et all (|2007d ). both these 
sources exhibit a strong spectral variation with spin phase. 
In the case of 4U 0142-1-614 the spectrum switches from be- 
ing very hard to very soft within a 0.1- wide phase interval 
and XMM-Newton data reveals a discontinuity, between 2 
keV and 3 keV, which can be interpreted as a curved com- 
ponent and is most appea rent within phase interval 0.7-0.9 
l|Den Hartog et aLll2008al '). In order to assess this scenario, 
a more detailed investigation of the pulse resolved spectra 
is necessary; this requires the introduction of the viewing 
geometry in our fits and will be presented in a forthcoming 
paper. 



4.2 Correlations 

We run a number of Spearman's rank correlation tests, in 
order to look for possible links among the observed proper- 
ties of the sources in our sample and the model parameters. 
In particular, we checked for correlations between the 1-10 
keV luminosity, I/i-io keV or the spin-down value of the mag- 
netic field, B, and each of the NTZ parameters, kT, Ptuik, 
and A<^. The values of the parameters are those obtained 
from the fit of the XMM-Newton data only (Tables [l] [H 
[S]). The source distances and the values of B have been taken 
from the compilation in iRea et al.l (|2008l ). The only signif- 
icant correlations which emerged are those between Pbuik 
and Li_iokeV and Ptuik and B. Both show a deviation from 
the null hypothesis probability of ~ 93%. Although the sig- 
nificance level is lower, 81%, a correlation between kT and 
ii-iokcv seems also to be present. The correlation between 
i'l-iokev and Pbuik, which is direct, can be explained tak- 
ing into account that an increase of Pbuik results in a larger 
energy gain of the photons per scattering and hence in a 
hardening of the spectrum, which translates into a higher 
luminosity. A similar argument applies to the Li_iokcV-fcr 
correlation, which is again direct. An increase of the sur- 
face temperature implies an increase of the flux of primary 
photons and again of the observed luminosity. The pbuik- 
B correlation m irrors that between Pt and B reported in 
iRea et al] l|2008l ). and, as discussed there, is not of immedi- 
ate interpretation. 

As discussed in iNobih. TuroUa fc Zan3 (12008^ 1 ■ the 
modulation of the X-ray flux may be due to the anisotropy 
of the magnetospheric charge density (which is lower along 
the poles), to a patchy surface temperature distribution, or 
to a combination of both effects. As mentioned in ^ in 
NTZ models the normalization constant (norm) divided by 
kT'^ is proportional to the emitting area on the star sur- 
face. To check if the observed modulation is related (also) 
to the presence of a hotter region on the surface, we looked 
for a correlation between the pulsed fraction and the ratio 
normx D'^/kT^, where D is the source distance. The pulsed 
fraction is defined as the semiamplitude of the sinusoidal 
function that best-fits the lightcurve, in the same energy 
range. We run the test on the entire sample of sources, ex- 
cluding the first observation of IE 1547.0-5408 , for which 
only an upper limit of the pulsed fraction is available. No 
significant correlation was found. However, when the sam- 
ple is restricted to AXPs only, a negative correlation (i.e. the 
pulsed fraction increases when the area decreases) emerges 
at the 91% confidence level. To verify this, we run again the 
test including the data from a set of five more recent ob- 
servations of IE 1547.0-5408 (Bernardini et al., in prepara- 
tion) and found that the correlation is still present, although 
the significance level slightly decreases to 74% (see Fig. IS}. 
Taken face value, this results points towards a localized emit- 
ting area in AXPs while the surface temperature would be 
more uniform in SGRs, where the modulation is produced 
mainly by scatterings in the magnetosphere. We note that 
the presence of a (time varying) hot spot has been reported 
in some transient AXPs, but the existence of a general pat- 
tern as the one discussed above needs a larger sample to be 
confirmed. 

Finally, we checked if there is a correlation between A(jf> 
and Pbuik- The motivation for this is that both parameters 
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are responsible for the formation of the high-energy tail. 
This may introduce a redundancy in the model parameters: 
in principle the fit might be not unique, since quite similar 
spectra may be obtained with different combinations of Ac/) 
and pbuik (e.g. low Ptuik and large or the opposite). We 
did not find any significant correlation between these two 
quantities. 

4.3 The broadband spectrum in the 1-100 keV 
range. 

As discussed in ^ a spectral decomposition of the kind 
NTZ+PL successfully reproduces the whole spectrum of 
IRXS J1708-4009 , IE 1841-045 , and SGR 1900-^14, up to 
200 keV. In this case, however, the best fit parameters of the 
NTZ model differ from those found from the fit of XMM- 
Newton data up to 10 keV. In particular, while the tem- 
perature of the seed blackbody is almost unchanged, the 
value of Pbuik is always considerably reduced (from 0.3-0.5 
down to 0.1-0.2) and the twist angle less constrained. This 
indicates that in the broadband fit most of the hardening is 
accounted for by the additional PL component that, in the 
case of IE 1841-045 and SGR 1900+14, starts to dominate 
the spectrum at energies as low as ~3-4 keV. Of course, 
this kind of spectral decomposition is certainly possible and 
may mimic a scenario in which the hard X-ray and soft X- 
ray emissions are due to two separate processes: for instance 
soft 7-rays may be produced in a twisted magnetosphere 
by thermal bremsstrahlung emission from the surface region 
heated by returning currents, or synchrotron emission from 
pairs created higher up (^100 k m) in the magnetosphere 
i|Thompson fc Belobodorwll2005l ). 

In principle, however, spectra produced by resonant cy- 
clotron upscattering of soft photons are expected to de- 
velop po wer law tails that exte n d up to much higher 
energi es l|Baring fc Harding! l2007l : iNobih. TuroUa fc Zand 
l2008al ). Therefore, it is well possible that our model can 
consistently explain the whole spectral energy distribution 
in the 1-200 keV range. If this is the case, the slope of the 
hard power law tail would be mainly dictated by the prop- 
erties of the magnetospheric electrons responsible for the 
upscattering, i.e. their density and velocity distribution. Ob- 
servationally the hard X-ray properties of AXPs and SGRs 
are quite different: while in the range 4-200 keV the spectra 
of SGRs are dominated by the non-terhmal component and, 
in the canonical model, well reproduced by a single unbro- 
ken PL, AXPs show a sort of turn over and have a non- 
thermal tail which is initially softer, up to ~ 5 — 8 keV, and 
then flattens at higher energies. In this picture, IE 1841-045 
seems an exception: its non-thermal emission is well repro- 
duced by a sing le PL and is therefore more SGR-like (see 
iRea et al] |2008| . and references therein). These properties 
suggest that, within the resonant Compton scattering sce- 
nario, AXPs spectra require different electron populations. 
On the other hand, a single electron distribution as for in- 
stance the one we used in the 3D computations presented 
here, may account for the broadband emission of SGRs or 
SGR-like sources. 

In order to test this possibility, we can not use directly 
our XSPEC table of models, since it has been computed 
assuming that magnetic scattering is conservative in the 
electron frame and, for self-consistency, spectra were trun- 



cated at 15 keV (see §[2] for a discussion). Instead, we re- 
computed the model that best-fits the 1-10 keV spectrum of 
SGR 1900-1-14 (which parameters are reported in Table [T|, 
this time by extending the computation to the whole 1- 
200 keV energy range and using a fully relativistic version of 
the Monte Carlo code w hich incorporates the complete QED 
scattering cross section l|Nobih. Turolla fc Zaneir2008bl ). We 
also repeated the computation by using slightly different sets 
of parameters, all within 3a from their best-fitting values 
(errors at la are reported in Table [ij . As expected, these 
new spectra differ from those computed in the Thomson 
limit only at very high energies (> 100 keV). Finally, we 
added to each curve a free normalization constant in the 
INTEGRAL range, to account for inter-calibration uncer- 
tainties between the two satellites. We then fitted the IN- 
TEGRAL data by leaving only the intercalibration constant 
as a free parameter. Results are reported in Fig.[7l where the 
two curves correspond to: i) the best fitting parameters re- 
ported in Tableland ii) to a different model in which (5buik, 
kT and Acp have been increased up to their upper limits 
within the 3a confidence range from the XMM-Newton fit. 
We found that in both cases the fit to INTEGRAL data 
is excellent (reduced ~ 1-14 and 1.07, respectively, for 
7 degrees of freedom), but the two model normalizations, 
in the XMM-Newton and INTEGRAL ranges, differ by a 
factor 8.8 and 2.8 respectively, too large to be attributed 
to intercalibrations uncertainties. Not surprisingly, similar 
large factors are also found when trying to fit the XMM- 
Newton and INTEGRAL data in the range 6-200 keV with 
a simple power law model. On the other hand, this simple 
test proves that the spectral slope of our model in the 20- 
200 keV range is the same as that of the INTEGRAL data. 
Similar considerations hold in the case of IE 1841-045 . Al- 
though a conclusive answer requires a direct fitting of the 
combined XMM-Newton and INTEGRAL data with a self- 
consistent Compton model, which is beyond the scope of 
this paper, we regard this preliminary finding as promising. 

4.4 Caveats and future developments 

We caveat that the models presented here are based on a 
number of simplifying assumptions. First, they are based 
on a globally twisted dipole model, that only gives an ideal- 
ized representation of the magnetic field topology. There are 
now both theoretical (Bclobo rodov 2009. ) and observational 
(a certain degree of hysteresis in th e long-term evolution 
of SGR 1806-20, IWoods et all l2007l : the long-term evolu- 
tion of the thermal component of t he transient AXP XTE 
J1810-197. IPerna fc Gotthelj l2008l : iBernardini eraLll2008D 
motivations in favour of a picture in which the twist affects 
only a limited portion the magnetosphere, typically the po- 
lar region. Furthermore, phase resolved spectroscopy of the 
two AXPs IRXS J1708- 4009 and 4U 0142+61 in the IN- 
TEGRAL energy range (jPen Hartog et al.|[2008al lbl) shows 
dramatic spectral changes with the spin phase. It has been 
recently suggested that a resonant scattering model in which 
the magnetic field is locally twisted ma y catch the essential 
features of this behaviour (|Pavan et al.i i2009l). 

A further point is the nature and velocity distribution 
of the magnetospheric charges. The NTZ model assumes 
the presence of mildly relativistic electrons, moving at con- 
stant velocity (which is a model parameter) along the closed 
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field lines. Currents flowing in t he maRnetosphere of a mag- 
netar have been investigated by iBeloborodov fc Thompson! 

who concluded that the magnetosphere is populated 
by pairs with a Lorentz factor 7 ~ 1000. However, more 
recent calculations (Nobili, TuroUa & Zane, in preparation) 
show that in the region where B > 2Bq Compton losses 
efficiently slow down electrons, limiting pair creation to a 
small region close to the star surface, and allowing for the 
presence of mildly relativistic particles along much of the 
circuit. 

Finally, angle of view effects have not been accounted 
for: the emerging spectrum is simply computed by inte- 
grating over the whole sky at infinity. This corresponds 
to the case in which the star is an aligned rotator, i.e 
the spin and magnet i c axes coincide. As discussed in 
iNobUi. TuroUa fc Zanj | 2008al ) . in order to treat the more 
general case in which the spin and magnetic axes are not 
aligned, we also produced a second XSPEC atable model 
by introducing two angles, x S'^d Ci which give, respectively, 
the inclination of the LOS and of the dipole axis with respect 
to the star spin axis. This allows us to take into account for 
the star rotation and hence derive pulse shapes and phase- 
resolved spectroscopy. The spectral model, ntzang.mod, has 
six free parameters {Pbuik, A^, log kT, x, C plus a normaliza- 
tion constant), i.e. two more than the one used here. Given 
that the NTZ model produces a very good fit for the phase- 
averaged spectra, the inclusion of two further parameters 
(i.e. the two angles) is not statistically required, as we tested, 
and will leave x ^-nd ^ unconstrained. On the other hand, 
having the possibility to infer the viewing geometry may 
be useful when fitting different outburst states in transient 
AXPs or when combining information that can be obtained 
by fitting simultaneously phase-resolved spectra, or indepen- 
dently from the study of the pulse profile. Further work in 
this direction is under way and will be presented in separate 
papers. 



5 CONCLUSIONS 

By considering a large sample of magnetars, we found that 
resonant Compton upscattering by a population of mildly 
relativistic electrons {Pbuik ~ 0.1 — —0.7) can reproduce 
the pulse averaged spectra in the range 1-10 keV. At var i- 
ance with the ID RCS model adopted in iRea et al. I (|2008l ). 
the approach used in the present investigation consistently 
accounts for the bulk motion of magnetospheric electrons 
which results in a more efficient comptonization of seed 
thermal photons. This has two main consequences: i) the 
required values of the optical depth are lower that those 
found using the ID RCS, and ii) NTZ spectra, being in- 
trinsically harder, successfully reproduce also the SGRs 
power-law tail below 10 keV. We found a significant cor- 
relation between the 1-10 keV source luminosity and both 
Pbuik and kT. This is indeed expected in the resonant cy- 
clotron scattering model and further supports this scenario 
to explain the high energy emission from magnetars. More- 
over, when restricting to AXPs only, we find hints for a 
negative correlation between the pulsed fraction and the 
emitting area. If confirmed, this suggests the presence of 
a strong thermal gradient on the star surface, which may 
also be responsible for the slightly different spectrum of 



4U 0142+614 and IE 2259-1-586 , the only two sources for 
which we could not find a satisfactory fit. Anisotropic sur- 
face thermal distributions may arise in the presence of large 
crustal magnetic fields, as expected in magnetars, because 
heat is preferentially transferred a long the field, resulting in 
small and confined hot caps (e.g. [Geppert. Kueker fc Page! 
200(i; IPons. MiraUes fc Geppertll2009l ). Moreover, the mag- 
netar surface is also expected to be further heated during 
bursting activity or by returning currents; in both cases the 
heat deposition can be substantially anisotropic. 

Finally, our conclusions regarding the modelling of the 
whole SED distribution up to ~ 200 keV are still com- 
patible with various scenarios. In the case of IE 1841-045 , 
IRXS J1708-4009 , and SGR 1900+14, a double component 
NTZ+PL model fits the combined XMM-Newton and IN- 
TEGRAL data, but it requires a NTZ spectrum quite soft 
and most of the emission dominated by the additional PL. 
This is compatible with a scenario in which soft X-ray and 
hard X-ray emission are ascribed to independent compo- 
nents. On the other hand, there is also the possibility that 
models of resonant upscattering can successfully describe 
the whole SED including the hard X-ray tail (eventually by 
invoking more than one electron population for sources with 
a spectral turnover at high energies). A conclusive answer 
requires a direct fitting of the combined XMM-Newton and 
INTEGRAL data, with all parameters left free and a new 
XSPEC table of models computed by using the QED cross 
sections. Further work on these topics is in preparation. 
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IRXS J1708-4009* 


IE 1841-045 


SGR 1900+14 


CXOU JOlOO-7211 


SGR 1627-41 


Parameters 




1.45±0.08 


2.5±0.1 


3.74±0.15 


< 0.05 


12±1 


kT (keV) 


0.47±0.09 


0.50±0.06 


0.45±0.04 


0.35 ±0.02 


0.64±0.11 


Pbulk 


0.34±0.04 


0.43±0.05 


0.46±0.05 


0.18±0.03 


0.7±0.1 


A(f> 


0.49±0.15 


0.47±0.04 


0.45±0.03 


1.7±0.8 


1.3±0.4 


NTZ norm 


0.35±0.04 


0.20±0.02 


0.08±0.02 


0.004±0.001 


0.005±0.001 


Flux (1-10 keV) (2.6 ±0.2) X 10"" 


(2.1 ±0.2) X 10"" 


(3.9 ±0.3) X 10-12 


(3.3 ± 0.2) X 10-13 


(3.3 ±0.2) X 10-13 


Xi (dof) 


0.97 (197) 


1.04 (152) 


0.99 (135) 


1.21 (101) 


1.16 (81) 



Table 1. Spectral Parameters of IRXS J1708-4009 , IE 1841-045 , SGR 1900+14 , CXOU JOlOO-7211 , and SGR 1627-41 obtained by fitting 
the XMM-Newton data with an NTZ model. The fit has been restricted to the 1-10 keV range, with the exception of CXOU JOlOO-7211 
and SGR 1627-41 for which we used the 0.1-10 keV and 2— lOkeV band, respectively. Errors are at la confidence level, reported fluxes 
are absorbed and in units of ergs-^cm-^. N/^ is in units of 10^^ cm-^ and we assumed solar abundances from Lodders (2003). See also 
Fig-IU Eland §[3] for details. 

*: source slightly variable in flux and spectrum, see text for details. 





IE 1547.0-5408 




IE 1048-5937 




Parameters 


2006 2007 


2003 


2005 


2007 




4.6±0.13 




0.66±0.02 




kT (keV) 


0.38±0.01 0.56±0.01 


0.60±0.06 


0.56±0.08 


0.69±0.02 


0bulk 


0.15±0.05 0.4±0.1 


0.14±0.03 


0.16±0.07 


0.13±0.04 




1.14±0.08 0.4±0.1 


1.9±0.2 


2.0±0.2 


1.9±0.1 


NTZ norm 


0.02±0.01 0.082±0.005 


0.10±0.03 


0.07±0.01 


0.21±0.01 


Flux (1-10 keV) (3.3 + 0.2) X lO-i^ (3.0 ± 0.2) x 10" 


■12 (1.1 ±0.2) x 10-11 ( 


8.3 ± 0.4) X 10-12 


(3.0 ±0.2) X 10-11 


Xl (dof) 


1.11 (164) 




1.22 (515) 





Table 2. Spectral Parameters of IE 1547.0-5408 and IE 1048-5937 in different emission states, obtained by fitting the XMM-Newton 
observations with an NTZ model. The fit has been restricted to the 1-10 keV range. Errors are at la confldence level, reported fluxes 
are absorbed and in units of ergs-icm-2. Nj:/ is in units of 10^2 cm— 2 and we assumed solar abundances from Lodders (2003). See also 
Fig.E] 



SGR 1806-20 


Parameters 


2003 


2004 


2005 






9.1±0.5 




kT (keV) 


0.69±0.10 


0.79±0.10 


0.74±0.15 


f^bulk 


0.48±0.05 


0.52±0.05 


0.51±0.05 


A<j) 


1.8±0.8 


1.7±0.7 


1.0±0.9 


NTZ norm 


0.11±0.05 


0.23±0.05 


0.12±0.08 


Flux (1-10 keV) (1.2 ±0.1) X 10" 


-11 (2.6 ±0.2) X 10" 


"11 (1.3 ± 0.2) X 10-11 


Xi (dof) 




0.98 (288) 





Table 3. Spectral Parameters of SGR 1806-20 in different emission states, obtained by fitting the XM M-Newton observations with an 
NTZ model. The three datasets have been taken before and after the giant flare of 2004 December 27 llHurlev et al.ll2005l : I Palmer et al.l 
l2005fl . with the intermediate observation taken about two months before the event. The fit has been restricted to the 1-10 keV range. 
Errors are at la confidence level, reported fluxes are absorbed and in units of ergs-lcm-2. Njj is in units of 1022 cm— 2 and we assumed 
solar abundances from Lodders (2003). See also Fig.|3] 
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IRXS J1708-4009* 


IE 1841-045 


SGR 1900+14 


Parameters 


NTZ+PL 


NTZ+PL 


NTZ+PL 


Nh 


1.67±0.03 


2.5±0.1 


3.9+0.2 


kT (keV) 

0bulk 

A4> 

NTZ norm 


0.39±0.04 
0.20±0.05 
1.90±0.04 
0.42±0.04 


0.50±0.06 
0.17±0.06 

> 0.72 
0.28±0.03 


0.43 ±0.06 
0.10±0.08 

> 0.38 
0.06±0.01 


r 

PL norm 


0.96±0.07 
(2.4 ±0.1) X 10--* 


1.46±0.08 
(2.1 ± 0.1) X 10-3 


1.87±0.09 
(1.1 ±0.1) X 10-3 


Flux (1-lOkeV) (2.6 ± 0.2) X 10"" 
Flux (1-200 kcV) (1.1 ±0.2) X 10""^ 


(2.1 ±0.2) X 10"" 
(1.1 ±0.2) X 10-1" 


(3.9 ± 0.3) X 10-12 
(1.5 ± 0.3) X 10-" 


Xi (dof) 


1.05 (205) 


1.12 (158) 


1.02 (140) 



Table 4. Spectral Parameters of IRXS J1708-4009 , IE 1841-045 , and SGR 1900+14 obtained by fitting the ~l-200keV XMM-Newton 
and INTEGRAL data with an NTZ+PL model. A constant function has been included to the modelling to take into account for 
intercalibration uncertainties (always < 10%). Errors are at Icr confidence level, reported fiuxes are absorbed and in units of ergs— ^cm- 2. 

is in units of lO'^^ cm— and we assumed solar abundances from Lodders (2003). See also Fig.|4]and §[3]for details. 
*: source slightly variable in fiux and spectrum, see text for details. 



IRXSJnOS-™ 1E1841-M5 SGR 1900+14 




i 10 1 2 5 lU 1 : j 10 

Energy (keV) Energy (keVj Eiicrgy(kcV) 



Figure 1. IRXS J1708-4009 , IE 1841-045 , and SGR 1900+14 : first row shows the spectra in Count/s/keV while in the second row we 
report the photon density plots of the modelling with the NTZ model. Only XMM-Newton data have been used, and the fitting has 
been restricted to the 1 — 10 keV range. See Table[T]and §[3]for details. 



X-ray spectra from magnetar candidates - III. 



11 





Figure 3. Same as in Fig.[T]for tiie AXPs IE 1547.0-5408 , IE 1048-5937 and for SGR 1806-20. See TablesE] |3] for details. 
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Figure 4. IRXS J1708-4009 , IE 1841-045 , and SGR 1900+14 : left column shows the spectra in Count/s/keV while in the right column 
we report the photon density plots of the modelling with the NTZ+PL model. Both XMM-Newton and INTEGRAL data have been 
used in the fitting. See Table|4]and §|3]for details. 




Figure 5. Spectrum from a single emitting zone on the star surface, computed as in Nobili, TuroUa & Zane (2008a, see Fig. 7 in that 
paper and text for details). Left: The emitting area is the equatorial strip ^ cos© ^ 0.25, 0^3"^ 7r/2. The LOS is at 0s = 90° 
and <I>s = 20° (dotted line), $s = 140° (dashed line) and <I>s = 220° (dash-dotted line). The solid line represents the seed blackbody. 
Right: The emitting area is the polar cap 0.75 ^ cos 8 ^ 1, ^ <J> ^ 27r. The LOS is at ©s = 19.5° (dash-triple dotted line), ©s = 63.5° 
(dashed line), ©s = 90° (dash-dotted line), ©s = 123.5° (dotted line) and ©j, = 160.5° (lower solid line). The upper solid line is the seed 
blackbody spectrum. 
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Figure 6. The pulsed fraction vs a measure of the surface emitting area for the AXPs in our sample; see text for details. Five more 
recent observations of IE 1547.0-5408 have been also included in this plot (Bernardini et al., in preparation). 
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Figure 7. Combined XMM-Newton and INTEGRAL spectrum of SGR 1900+14 . The dashed line is the NTZ spectrum computed for 
the same values of parameters that best-fit the XMM-Newton data only (see Table [TJ, but using the QED cross sections in the Monte 
Carlo code. This spectrum is than used to fit the INTEGRAL data alone by leaving free only the value of the normalization. This gives 
= 1.14 for 7 d.o.f. and the ratio between the INTEGRAL and the XMM-Newton normalizations is 8.8. Solid line: same as before but 
with fShulk = 0.61, kT = 0.57 keV and A(f> = 0.54; all parameters are at 3cr from their best fitting values of the XMM-Newton spectrum. 
Here the reduced = 1-07 for 7 d.o.f. and the ratio between the INTEGRAL and the XMM-Newton normalizations is 2.8. 



